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INTRODUCTION  AND  SUMMARY 


An  objective  of  the  CIRRIS  Program*  is  to  obtain  measurements  that 

(a)  define  the  Earth's  limb  as  an  IR  spectral  background,  and 

(b)  provide  Information  on  photochemical  processes  in  both  the 
quiescent  and  disturbed  atmosphere. 

Interpretation  of  the  limb  data  In  terms  of  atmospheric  photochemistry  will 
include  mathematical  inversions  of  the  data,  i.e,,  transformation  of  the 
measured  limb  radiance  profiles  into  vertical  distributions  of  concentration 
and  temperature. 

CIRRIS  will  obtain  spectrally  resolved  Earth  limb  radiarce  profiles  up 
to  several  hundred  kilometers  tangent  height.  Quantities  that  can  be 
inferred  for  the  nonequilibrium  atmosphere  (H  'S'  70  km)  through  an  inversion 
of  the  data  include  ground  state  and  excited  state  species  concentration  and 
rotational  temperature.  However,  limb  inversion  techniques  applicable  to 
spectrally  resolved  data  and  to  the  non-LTE  altitude  regime  have  not  been 
studied  extensively.  These  differ  significantly  from  limb  inversion  methods 
developed  for  the  lower  atmosphere  (e.g.,  Gille  and  House,  1971). 

The  purpose  of  the  present  study  was  to  develop  an  algorithm  applicable 
to  a  special  nonequilibrium  case,  namely,  the  recovery  of  the  distribution 
of  excited  nitric  oxide  and  kinetic  temperature  above  100  kilometers  from 
limb  spectral  radiances  measured  in  the  fundamental  (Av  »  1)  band  at  5-3  pm, 
This  is  a  special  case  because  ground  state  NO  is  optically  thin  for  tangent 
paths  above  ~ 1 00  km  for  the  fundamental  band.  The  study  results  provide  a 
starting  point  for  development  of  algorithms  applicable  to  the  non-optical ly- 
thln,  nonequilibrium  case,  which  will  be  required  in  analyzing  CIRRIS  data. 

The  developed  algorithm  was  applied  to  NO  spectral  radiance  profiles 
measured  by  SPIRE  (Spectral  Infrared  Rocket  Experiment).  The  excited  NO 


*CIRRIS  Is  an  acronym  for  Cryogenic  J_nf raRed  Radiance  Jjistrumentation  for 
Shuttle.  The  CIRRIS  Program  is  descVibecT  by“Ahmadj  ian,  et  at  (1980). 


densities  were  used  to  compute  atmospheric  radiative  cooling  rates  for  the 
5.3  pm  band,  and  the  de-excitation  rate  constant  for  the  reaction 
NO*  +  0  ■+•  NO  +  0,  These  results  are  discussed  in  Sections  3  and  k.  Section  2 
describes  the  inversion  algorithms.  The  SPIRE  instrument  and  data  it  obtained 
are  described  by  Nadile  et  al  (1977)  and  Stair  et  at  ( 1 981 ) . 

The  solution  for  NO  excited  density  and  the  computed  cooling  rates  and 
de-excitation  rate  constant  are  in  good  agreement  with  other  published 
measurements  and  models.  The  solution  for  temperature  is  limited  by  noise  in 
the  SPIRE  spectral  data. 


SECTION  2 

INVERSION  ALGORITHMS 

The  geometry  of  the  limb  observations  Is  shown  In  Fig.  la.  The  band 
radiance  at  the  sensor  due  to  excited  NO  molecules  in  path  segment  dz  Is 

dN(HT)  -  (AAtt)  d2  n*  (z)  ,  (1) 

where  HT  Is  the  tangent  height  of  the  observation,  n*(z)  is  the  density  of 
the  excited  molecules  at  altitude  H(Hy,z)  and  A  Is  the  Einstein  A-coeff iclent 
for  the  band  (NO  fundamental).  We  have  assumed  that  reabsorption  by  ground 
state  NO  over  the  tangent  path  can  be  neglected.  The  band  radiance  due  to  the 
entire  path  can  be  expressed  as  an  integral  over  altitude: 

2A  F  * 

"(ht>.  -  $  J dH  -sir- n  (H)  •  (2) 

ht 

A  form  suitable  for  numerical  evaluation  or  inversion  is 

X 

N(Hj)  s  Nj  ■  (2A/Att)  ^  A,_m  n*  ,  (3a) 

1  -J 

where  the  A  j  j  are  distances  (Az's)  defined  in  Fig.  lb.  Given  a  limb  profile 
of  the  band  radiance  (set  of  measurements  Nj) ,  one  can  invert  Eq.  (3a)  to  find 
a  set  of  values  n*  that  approximates  n*(H).  The  set  Nj  is  obtained  by  Inte¬ 
grating  the  SPIRE  limb  radiance  spectral  data  over  the  NO  fundamental  band. 

Equation  (3a)  will  give  a  good  approximation  to  the  observed  limb 
radiance  profile  If  n*  does  not  vary  appreciably  over  the  vertical  distance 
AH  resolved  by  the  sensor;  note  that  AH  determines  the  weights  Ajj.  If  the 
variation  Is  appreciable,  a  smaller  AH  should  be  used  in  evaluating  Eq.  (3a), 
and  the  result  should  be  convolved  with  a  function  representing  the  sensor's 
spatial  response  In  the  verttcal. 

The  limb  profile  of  spectral  radiance  measured  by  a  device  like  the 
SPIRE  CVF  spectrometer  can  be  expressed  In  the  form 
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I 

"<Vk>  5  NJk  '  "u  •  <3"> 

1  -J 

where 

"ill  ■  nTRik  <« 

and 

R[T(Hj),Xk]  s  RJk  -  2Si[T(H,)]  C(At,Xk)  ,  (5) 

l 


where  X^  denotes  the  centers  of  the  rotational  lines  comprising  the  band,  S £ 
are  the  line  Intensities  (determined  by  the  NO  rotational  temperature  at 
altitude  Hj),  and  C ( , Xk)  Is  the  spectral  response  of  the  Instrument  at 
wavelength  Xk  due  to  a  line  (unit  impulse)  at  X  «  X£.  The  response  function 
C  and  line  Intensities  S£  are  normalized  such  that 


Thus,  Eqs.  (3b),  (k)  and  (5)  merely  replace  the  band  source  function  A  of 
Eq.  (3a)  by  the  spectral  source  function  AR,  which  depends  on  rotational 
temperature.  Evaluation  of  the  normalized  apparent  spectral  emission  R  is 
simplified  by  the  fact  that  the  Doppler  broadened  lines  are  very  narrow 
compared  to  the  response  function  C;  they  can  be  treated  as  delta  functions. 
However,  the  instrument  response  function  C  (which  is  dominated  by  the  spectral 
response  of  the  CVF)  depends  on  both  X^  and  Xk,  not  on  Xk~X£;  i.e.,  R  cannot 
be  represented  as  the  convolution  of  a  single  function  with  a  series  of  impulses. 


Clven  a  data  set  N,.  versus  j  for  fixed  k  (limb  profile  of  spectral 


radiance  at  fixed  spectrometer  setting  X.),  one  can  invert  Eq.  (3b)  in  the 

ft  ft  " 

same  manner  as  Eq.  (3a)  to  obtain  n.  versus  i.  Then  R..  can  be  found  from 

^  I K  I  ^ 

Eq.  (A),  using  the  solution  nj  obtained  by  inverting  the  band  radiance 


profile.  By  repeating  the  process  at  other  wavelengths  one  obtains  Rjk 


versus  k  and  i,  I.e.,  the  normalized  spectral  distribution  R  versus  X 


for  different  altitudes.  Finally,  each  of  these  recovered  spectral  distribu¬ 
tions  can  be  compared  to  ones  computed  from  Eq.  (5)  for  various  rotational 
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temperatures  T,  to  determine  the  temperature  that  gives  the  best  agreement, 

The  result  is  a  solution  for  the  vertical  distribution  of  the  NO  rotational 
temperature  (which  is  the  same  as  the  kinetic  temperature,  except  at  alti¬ 
tudes  greater  than  —  1 60  km). 

Note  that  this  procedure  for  recovering  ro.ational  temperature  is  not 
restricted  to  the  case  of  low  spectral  resolution.  However,  if  the  sensor 
can  resolve  individual  lines,  It  may  be  desirable  to  use  integrated  line 
radiances  rather  than  spectral  radiances,  and  to  redefine  N..  and  R.. 

A  JK  IK 

accordingly.  The  SPIRE  data  used  to  recover  the  temperature  profile  con¬ 
sisted  of  10  spectrum  values  with  wavelength  spacing  between  0.0$  ym  and  0.  urn. 

Equations  (3a)  and  (3b)  were  inverted  using  a  so-called  direct  nonlinear 
technique.  The  method  starts  with  an  initial  guess  for  the  solution  n,  or 
n.^,  and  uses  the  equations  to  compute  the  corresponding  limb  profile  Nj  or 
Nj^.  Ratios  between  the  measured  and  computed  profile  values  are  used  to 
relax  the  initial  guess.  The  process  is  iterated  until  the  rms  difference 
between  computed  and  measured  limb  radiances  is  smaller  than  some  preassigned 
value,  which  might  be  the  system  noise.  Solutions  obtained  by  such  a  method 
are  generally  non-unique,  and  the  number  of  values  comprising  the  solution 
may  exceed  the  number  of  observations.  The  efficiency  of  the  inversion  proce¬ 
dure  and,  to  some  extent,  the  quality  of  the  'nlutlon,  depends  on  the 
particular  relaxation  method  employed. 

2.1  Relaxation  Method 


Equation  sets  (3a)  and  (3b)  evaluated  in  the  nth  iteration  have  the  form 


m  y  B<«>  +  y  >> 
WUSl  W21*2 

+  W  s(n) 

31s3 

+  . . . 

W  s(n) 
W 1 2S2 

+  W  s(n) 
**22s  3 

+  ... 

(7) 

(-.) 

3 

K 

W  s(n) 
Wl  3S3 

+  ... 

The  choice  of  spectral  resolution  for  a  limb  inversion  sensor  involves  many 
issues,  including  available  signal-to-noise  and  the  possibility  of  spectral 
interferences  from  other  atmospheric  species.  Low  to  moderate  resolution 
may  be  optimum  in  some  circumstances. 
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where  the  s;  1  represent  the  Iterated  guess  for  the  solution  profile,  W 

,(n)  mj 


’1 


ere  the  weights  (2A/4ir)A  , ,  end  0,  '  ere  the  corresponding  computed  limb 

J  /Q\ 

redtences  or  spectre!  radiances.  Let  0j  denote  the  limb  observations. 
The  relexetion  procedure  used  in  the  present  study  can  be  described  as  a 
seeling  end  averaging  operation: 

a.)  For  the  Jth  tangent  height  (jth  row  of  Eg.  (7)),  calculate  the 
quantity  (oj°^/Ojn^ ) *.  Multiply  each  s^  of  the  Jth  row  by 
thi s  quantity. 


b.)  After  scaling  each  row  In  this  manner,  obtain  s 


(n+1 ) 
l 


as  the  sum 


of  the  elements  of  column  1  divided  by  the  sum  of  the  weights 

.(n+1) 

'  I 


The  s,  ",/  are  the  new  solution,  or  final  solution 


i n  co 1 umn  1 . 

If  the  convergence  criterion  Is  satisfied. 

This  Is  the  scaling/averaging  operation  proposed  by  Halem  and  Chow  (1976) 
for  the  nadir  Inversion  problem,  except  that  in  the  present  case  the  informa¬ 
tion  channels  correspond  to  tangent  height  rather  than  spectral  subbands. 

The  formal  relaxation  equation  for  the  limb  geometry  is 


(n+1)  „  (n) 


\' 

j-i 


‘i-j+ij  ; 


(8) 


The  optimum  value  for  the  exponent  K  is  determined  by  experimentation.  The 
value  Koj  1.5  was  found  to  yield  the  maximum  rate  of  convergence  In  the  inver¬ 
sions  of  SPIRE  NO  data. 


The  measured  limb  radiance  profile  has  upper  and  lower  bounds  on  tangent 
height  corresponding  to  subscripts  j  ■  1  and  J,  respectively.  The  upper 
bound  is  determined  by  the  instrument  sensitivity  (noise)  or  the  upper  limit 
of  the  spatial  scan.  It  Is  for  this  reason  that  the  limit  i '  in  the  summa¬ 
tions  in  Eq.  (8)  is  equal  to  the  smaller  of  the  Indices  i  and  J.  Note  that 
the  number  of  solution  values  I  can  be  greater  than  the  number  of 
observations  J, 
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The  weights  Amj  in  Eqs.  (3a).  (3b)  end  (8)  are  given  by  (see  Fig.  lb) 


A  .  *  z  .  -  z  ... 

ipj  nj  m- 1 ,  j 


(9) 


where 


20j  =  °  ' 


1/2 


(10) 


zmj  *  {[R  +  Hj  +  (m  -  1  -H/2)AH]2  -  (R  +  H.)2}  ;  m-1,2 . 


and  R  is  the  Earth's  radius.  Note  that  z^j  can  be  approximated  by 

zmj  a  [AN(R  +  HT)31/2(2m-  1),/2  =  zm  (11) 

where  ifT  Is  an  average  tangent  height.  The  maximum  relative  error  In 
due  to  this  approximation  is  about  one  percent  over  the  range  of  altitudes 
involved  in  the  NO  inversions.  Thus,  the  weights  Afflj  which  appear  in  both 
numerator  and  denominator  of  the  relaxation  formula,  Eq.  (8),  can  be  replaced 
by  the  much  simpler  quantity  w^: 

Wj  ■  1 

w  -  (2m-  l)l/2  -  Um-3)’/2!  m-2,3 .  (12) 

m 


This  approximation  was  used  for  relaxing  the  guess  solution.  The  exact 
equations  (S)  end  (10)  were  used  in  evaluating  the  limb  radiances  for  the 
iterated  guess. 


2.2  Analytical  Solutions  for  a  Special  Case 

Studies  dealing  with  the  reaction  chemistry  and  thermal  budget  of  the 
upper  atmoH-nere,  such  as  Gordiets,  et  at  ( 1 980)  and  Kockarts  (l98o)»  have 
derived  or  made  use  of  models  in  which  the  density  of  excited  NO  varies  nearly 
exponentially  with  altitude  above  H  a  150  or  200  km.  It  can  be  shown  that 
the  function  N(H.,)  is  approximately  proportional  to  n*(H)  when  the  latter  is 
exponential;  i.e.,  the  limb  radiance  profile  is  approximately  an  exponential 
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b.k- 


with  the  same  decay  constant.  More  generally,  the  contribution  of  an 
arbitrary  atmospheric  layer  to  the  limb  radiance  at  some  tangent 

height  £  Hj  can  be  expressed  in  an  approximate  closed  form;  that  is,  in 
terms  of  an  exponential  and  the  error  function.  Tne  analytical  solutions 
are  given  in  Appendix  A. 

The  SPIRE  integrated  limb  spectra  are  roughly  exponential  at  the  higher 
limb  tangent  heights  (Hj  £  175  km)  if  the  noise  is  averaged;  the  noise 
becomes  significant  above  175  km.  SPIRE  data  for  the  high  altitudes  and  any 
other  measurements  that  exhibit  the  exponential  behavior  can  be  "instantly 
Inverted"  using  the  analytical  solution.  The  more  general  solution  given  in 
Appendix  A  is  useful  in  establishing  practical  limits  for  the  integrations 
required  in  the  inversions. 

2.3  Inversion  using  Fourier  Transforms 

The  geometric  approximation  used  to  simplify  the  relaxation  equation  and 
to  obtain  the  analytical  solutions  for  an  exponential  limb  radiance  profile 
can  also  be  used  to  show  that  Eqs,  (3a)  and  (3b)  can  be  expressed  approximately 
as  convolutions.  This  means  that  the  inversion  (for  the  optically  thin  case) 
can,  in  principle,  be  performed  as  a  deconvolution,  that  is,  by  Fourier  trans¬ 
form,  multiplication  and  inverse  Fourier  transform  operations.  The  method  is 
described  in  Appendix  B.  It  is  possible  that  the  method,  implemented  via  the 
FFT  algorithm,  is  more  efficient  than  the  direct  nonlinear  method  described 
above.  However,  we  have  not  tested  the  method.  Unlike  the  direct  nonlinear 
method  it  cannot  be  adapted  to  the  non-optical ly-thin  case. 
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SECTION  3 

LIMB  RADIANCE  DATA  AND  INVERSION  RESULTS 


!t  was  determined  that  a  reasonable  altitude  increment  AH  for  the 
numerical  integrations  (see  Fig.  lb)  is  2,5  km.  The  upper  part  of  Fig.  2 
shows  an  approximation  to  the  weighting  function  3z/3H  and  the  corresponding 
numerical  weights  evaluated  for  AH  -  2.5  km.  The  weighting  function  in  the 
figure  is  3z/3H  »  const,  x  (H-H.J.)  (see  Appendix  A),  and  the  numerical 
weights  are  the  integral  of  this  approximation  over  the  intervals  (Ht,  Hj  + 
0.5  AH),  (Ht  +  0.5  AH,  Ht  +  l.s  AH),  (Ht  +  1.5  AH,  H?  +  2.5  AH),  ...;  i.e., 

the  latter  are  given  by  Eq.  (12).  The  two  dashed  curves  show  the  maximum 
relative  variation  of  excited  NO  density  with  altitude  in  a  model  used  by 
Kockarts  (1980).  It  is  seen  that  an  integration  mesh  size  of  2.5  km  is 
adequate  for  computing  the  limb  profile  due  to  excited  NO,  if  the  actual 
gradients  are  no  larger  than  those  of  Kockarts'  model.  The  value  used  for 
A  =  „  ■  Aq^j,  the  thermally  averaged  Einstein  emission  coefficient 

for  the  v-  «*  1  transition  of  the  NO  fundamental,  is  10.78  sec  after 
Billingsley  (1975). 

The  SPIRE  instrument  resolved  appropriately  5  km  at  the  tangent  height, 
for  tangent  heights  between  120  and  160  km.  Therefore,  it  was  necessary  to 
degrade  the  spatial  resolution  of  the  limb  profiles  computed  by  Eqs.  (3a)  and 
(3b)  before  comparing  them  to  the  SPIRE  measured  profiles.  The  procedure 
consisted  of  convolving  the  Nj  with  the  three-point  numerical  comb  (1/4,  1/2, 
1/4),  which  represents  the  sensor  1 FOV  as  a  rectangular  (boxcar)  Function  of 
width  5  km.  The  effective  numerical  weighting  function  for  the  SPIRE  observa 
tions  is  shown  in  the  lower  part  of  Fig.  2.  This  was  obtained  by  convolving 
the  weights  in  the  upper  part  of  the  figure  with  the  numerical  IFOV  comb. 

Figure  3  shows  several  band  radiance  limb  profiles  obtained  by  inte¬ 
grating  the  SPIRE  spectral  t  ta  over  the  NO  5.3  pro  band.  Scans  4  and  9 
observed  the  limb  in  daylight.  In  scan  8,  tangent  heights  below  ~13Q  km 
are  In  darkness;  the  relatively  high  radiance  levels  between  90  and  120  km 
suggest  that  this  scan  may  have  observed  a  weak  aurora  (see  Stair  ei  at  1 98 1 ) 
The  spacing  of  the  data  in  tangent  height  is  approximately  four  to  eight  km. 
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WEIGHTING  FUNCTION  OR  NUMERICAL  WEIGHT 
(ARBITRARY  UNITS) 


Figure  2 


Approximate  weighting  function  for  optical ly-thin  limb  inversion 
and  the  corresponding  numerical  weights  for  the  Integration  mesh 
size  AH  «  2.5  km. 
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figure  3.  Umb  radiance  profiles  obtained  by  SPIRE  for  the 
NG  fundamental  band  (5. 08-5. hk  yn,'  . 
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The  SPIRE  noise  equivalent  band  radiance  is  approximately  2  x  10  °Wcrn  sr’ 
in  the  NO  band.  The  noise  is  evident  at  the  higher  tangent  heights. 

The  continuous  curve  fn  Fig.  4  Is  the  spectral  radiance  measured  by 
SPIRE  at  tangent  hefont  Hy  *  122  km  during  scan  9>  The  nonzero  baseline  is 
an  artifact  caused  by  saturation  of  the  detector  at  an  earlier  time  in  the 
CVF  spectral  scan  sequence.  The  baseline  effect  is  greatest  at  the  lower 
tangent  heights,  and  becomes  negligible  when  Hy  ?  130  km.  The  baseline  was 
assumed  to  be  linear  with  wavelength  and  was  estimated  from  spectral  radiances 
near  the  band  edges;  it  was  removed  before  the  data  was  inverted.  Data  for 
Hy  <  100  km  was  not  used  because  it  is  not  clear  that  the  correction  method 
is  valid  when  the  baseline  spectral  radiance  is  large.  Figure  3  shows 
uncorrected  profiles.  The  dashed  curve  in  Fig.  4  is  described  in  a  later 
section. 

3 . «  Results  for  n*(H) 

The  SPIRE  data  for  scan  9  was  Inverted  to  obtain  the  vertical  distribu¬ 
tions  of  excited  NO  density  n*(H)  and  kinetic  temperature  T(H).  The  measured 
profile  of  band  radiance  is  shown  in  Fig  5.  The  points  in  this  figure  were 
obtained  by  removing  the  nonzero  baseline  and  Integrating  the  resulting 
spectra  between  X  ■  5.0838  pm  and  5.6370  pm.  The  noisy  data  above  Hy  =  192  km 
were  replaced  by  the  indicated  exponential  profile.  The  corrected  and 
smoothed  data  were  interpolated  onto  a  uniform  grid  with  intervals  AHy  *  2.5  km 
between  Hy  *  100  and  250  km.  The  Integration  mesh  used  to  compute  iimb 
radiances  extended  to  400  km. 

Two  inversion  results  n*(H)  are  also  shown  In  Fig.  5.  Curve  (a),  which 
exhibits  excessive  vertical  structure,  was  obtained  by  iterating  until  the 
rms  relative  difference  between  the  observed  and  computed  radiances  was 
approximately  one  percent.  This  solution  required  34  iterations.  The  other 
solution,  continuous  curve  (b) ,  represents  an  rms  relative  difference  of 
4.1  percent,  and  required  six  Iterations.  The  result  for  five  iterations 
(not  shown)  had  a  relative  rms  difference  of  6.8  percent,  and  if  plotted 
would  be  barely  distinguishable  from  curve  (b) .  The  initial  guess  in  each 
case  was  n*(H)  ■  const.  *  lO^  cm  ^ . 
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Figure  5.  Corrected  and  smoothed  NO  limb  radiance  profile,  and  excited 

NO  density  profiles  n*(H)  obtained  by  Inversion,  The  smoother 
solution  (curve  b)  corresponds  to  a  four-percent  relative  rms 
difference  between  observed  and  computed  limb  radiance  profiles. 


The  sharp  structure  in  solution  (a)  near  175  km  Is  due  to  the  sharp 
feature  in  the  limb  radiance  profile  near  this  altitude,  which  is  probably 
a  noise  "spike".  Note  that  the  sensor,  as  it  scans  downward  in  tangent 
height,  tends  to  produce  a  radiance  profile  that  is  roughly  the  Integral  of 
the  n*(H)  profile;  that  is,  sharp  features  are  integrated  because  of  the 
sharp  lower  cutoff  of  the  weighting  function  (Fig.  2).  Thus,  the  solution 
can  display  unrealistic  S-shaped  features  that  resemble  the  derivative  of 
noise  spikes  in  the  measured  radiance  profile.  Solution  (b)  shows  that  noise 
has  a  smaller  effect  when  the  convergence  criterion  is  set  at  a  value  that  is 
closer  to  the  expected  rms  relative  noise  in  the  measurements.  For  the 
purpose  of  deriving  reaction  rate  constants  and  atmospheric  cooling  rates 
(Section  A),  we  used  solution  (b) ,  but  smoothed  this  result  in  the  region  of 
170  km.  The  dashed  curve  in  Fig.  5  represents  the  smoothed  portion  of  the 
solution. 

We  found  that  the  choice  for  the  initial  guess  affected  the  rate  of 
convergence,  but  had  virtually  no  effect  on  the  final  solution. 

3.2  Results  for  T(H) 

Careful  examination  of  the  SPIRE  spuctral  radiances  and  their  variation 
with  tangent  height  indicated  that  there  were  unknown  shifts  in  the  wave¬ 
length  scale  of  the  data.  The  shifts,  which  were  in  the  range  ±  0.03  pm, 
appeared  random  with  tangent  height  for  some  scans.  We  elected  to  use  the 
SPIRE  scan  9  data  for  the  inversions  because  it  had  an  apparent  shift  that 
was  a  relatively  slow  function  of  tangent  height.  The  wavelength  shifts  were 
not  large  enough  to  seriously  affect  the  band  radiances  calculated  from  the 
data . 


The  shifts  were  removed  by  the  following  procedure.  A  reference  tempera¬ 
ture  profile,  the  solution  n*(H)  obtained  from  the  scan  9  band  radiances,  and 
Eqs .  (3b),  14)  and  (5)  were  used  to  compute  a  limb  profile  of  spectral 
radiance  N(Hy,x).  The  computed  spectrum  for  *  122  km  (added  to  the  base¬ 
line  correction)  is  shown  by  the  dashed  curve  in  Fig,  4.  The  locations  of 
the  P  and  P  branch  maxima  In  the  computed  spectrum  relative  to  those  in  the 
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measured  spectrum  were  used  to  obtain  the  wavelength  shift  correction  to  the 
latter.  We  used  the  temperature  profile  given  by  the  Jacchia  (1977)  model 
for  an  exospheric  temperature  1050K  (based  on  the  U.S.  Standard  Atmosphere 
1976)  as  the  reference  temperature  profile. 

After  each  of  the  measured  spectra  were  shifted  in  this  manner,  a 
vertical  temperature  profile  was  recovered  using  the  rotational  temperature 
inversion  procedure  outlined  in  Section  2.  It  can  be  argued  that  the  solu¬ 
tion  has  been  constrained  by  the  necessity  to  assume  a  temperature  profile 
in  the  correction  of  the  unknown  wavelength  shifts.  There  is  surely  some 
measure  of  validity  in  this  argument.  However,  noise  in  the  spectral  data 
probably  had  a  much  larger  effect  in  the  solution  than  the  method  used  to 
correct  for  the  shifts,  and  would  probably  have  dominated  even  if  no  correc¬ 
tions  had  been  made. 

Figure  6  shows  the  SPIRE  scan  9  measured  spectral  radiances  for  tangent 
heights  of  148,  160  and  177  km.  Because  of  the  rapid  deterioration  of 
signal-to-noise  above  HT  ■  150  km  we  made  no  attempt  to  recover  the  tempera¬ 
ture  above  150  km  altitude.  Also,  the  data  used  for  the  inversions  was 
spectrally  smoothed  up  to  HT  ■  166  km,  and  was  discarded  above  this  tangent 
height.  We  assumed  that  the  spectral  distribution  of  observed  limb  radiance 
remains  constant  for  HT  £  166,  and  that  the  spectral  radiances  are  propor¬ 
tional  to  the  band  radiance.  This  assumption  implies  constant  temperature 
for  H  £  166,  but  should  not  have  a  large  effect  on  the  retrieved  temperature 
profile  below  H  ■  150  km. 

Some  of  the  solutions  RfT(Hj),Xk]  obtained  by  inversion  of  the  SPIRE 
data  are  shown  in  Fig.  7a.  The  ten  selected  wavelengths  Xk  are  5.1.  5.15, 
5-2,  5.25,  5-3,  5.35,  5-1*,  5-5,  5.6  and  5.7  ym.  Figure  7b  shows  theoretical 
R(T,\)  values  computed  from  Eq.  (5)  for  various  temperatures  in  the  range 
300K  to  700K.  The  generally  poor  agreement  between  the  set  of  recovered  R's 
and  computed  R's  is  a  result  of  the  poor  S/N  of  the  spectral  data,  and 
possible  wavelength  scale  "stretching"  in  the  data. 
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Figure  7.  Some  of  the  recovered  spectral  distributions  R[T(Hj),X^],  and 
theoretical  spectra  R(T,Xj<).  R  Is  the  normalized  apparent 
spectral  emission  per  unit  volume  at  rotational  temperature  T 
(see  Eq.  5).  and  has  units  ynfl. 


Figure  8  shows  the  recovered  temperature  profile  obtained  by  a  systematic 
comparison  of  the  solutions  R-k  and  a  large  library  of  computed  R(T,X) 
spectra.  Plotted  in  the  same  figure  is  the  temperatui  a  piovilw  given  by  the 
Jacchia  model.  The  solution  T(h),  although  extremely  noisy,  Is  in  rough 
agreement  with  the  Jacchia  profile  above  H  c;  120  km.  The  S-shaptd  temperature 
excursions  resemble  those  produced  In  the.  n’'(ll)  solut5on  by  noise. 
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The  temperature  profile  T(H)  obtained  by  Inversion  of  the  SPIRE 
NO  spectral  data,  and  a  profile  given  by  the  Jacchia  (1977) 


SECTION  It 

INTERPRETATION  OF  RESULTS 


KocKarts  (1980)  pointed  out  that  emission  by  the  NO  fundamental  band 
is  a  major  cause  of  heat  loss  In  the  upper  atmosphere.  He  gives  a  simple 
formulation  of  the  process,  wnicn  results  in  a  cooling  race  L  for  the 
5.3  ym  band  given  by 

L(H)  -  Ahcvn*(H)  (13) 

where  h  is  Planck's  constant,  c  Is  the  velocity  of  light  and  v  Is  the 
wavenumber  of  the  NO  emi ssion .  (Ahcv  st  4.1  x  !0  *  If  L  has  units  erg  cm  Js  1 
and  nA  has  units  cm~^) ,  The  solution  n*(H)  obtained  by  inversion  of  the 
SPIRE  data  was  used  to  calculate  L.  The  result  is  compared  in  Fig.  9  to 
cooling  rates  given  by  Kockarts  and  by  Gordiets  at  al  (1980).  The  three 
results  agree  within  a  factor  of  four  at  H  *  120  km,  and  within  a  factor  of 
two  above  140  km.  These  differences  are  consistent  with  expected  variations 
In  n  and  n*. 


A  steady  state  condition  on  the  four  important  processes  involving  NO 
in  the  upper  atmosphere  (spontaneous  emission  by  e-clted  NO,  excitation  and 
deactivation  of  NO  by  collisions  with  u  atoms,  and  excitation  of  NO  through 
absorption  of  Earthshine)  corresponds  to  the  concentration  ratio 


* 


+  R 

n  y  +  A 


(14) 


where 


n  «  the  concentration  of  ground  state  NO, 

K  *  exp(-hcv/kT)  a  exp ( “2698. 5°K/T) , 

Y  -  the  rate  of  excitation  of  NO  due  to  0  atoms, 

A  *  the  Einstein  emission  coefficient  ■  10.78  s  *, 

R  *  the  rate  of  production  of  excited  NO  due  to  absorption  of  Earthshine; 

R  *  irB  S  where 
v 

B  (T-)  -  the  Planck  spectral  emission  rate  corresponding  to  the 
v  "  effective  Earthshine  temperature  T^,  and 
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S  -  the  band  strength  of  the  5.3  ym  band  cs  503  x  10  crn/molec 
at  296K  (Rothman  et  al  1981). 

The  excited  concentration  n*  is  known  from  inversion  of  the  SPIRE  NO  data, 

If  models  are  adopted  for  n(H)  and  T (H ) ,  then  we  know  i|)(H)  and  K(H) ,  and  can 
solve  Eq.  (14)  for  y(H)  : 

y(H)  »  (A<j,-R)/(<|>-K)  (15) 

Over  the  altitude  range  for  which  Eq.  (15)  will  be  evaluated,  R  can  be 
neglected  in  comparison  to  A \l>.  The  de-excitation  rate  constant  kg  for  the 
reaction 

k0 

NO*  +0  +  H0+0  (16) 

can  be  obtained  from  y(H)  -  kg[0]  if  a  model  is  also  adopted  for  the  concentra¬ 
tion  of  oxygen  atoms  [0(H)]. 

Table  1  shows  the  results  of  the  calculation  of  k^.  The  kinetic  tempera¬ 
ture  profile  is  the  one  given  by  the  Jacchia  (1977)  model,  as  described 
earlier.  The  distribution  n(H)  =  [NO (H) ]  is  the  "median"  model  suggested  by 
Caledonia  et  al  (1980),  and  the  profile  for  [0]  is  taken  from  Jacchia  (1977). 

The  average  of  k_  over  the  altitude  range  1 30- 1 60  km  is  ~7.8  x  10  ' '  cm  3  s 
u  - 1 1 

which  compares  favorably  with  the  value  6.5  x  10  measured  by  Fernando  and 
Smith  (1979)  at  a  temperature  of  300K.  Glazer  and  Troe  (1975)  measured  the 
value  3-6  x  10  ^  cm  * s  '  at  2700K.  If  we  yse  the  "auroral"  model  for  [NO] 
given  by  Caledonia  et  al  (1980),  which  has  NO  concentrations  approximately 
five  times  higher  than  the  "median"  model,  we  obtain  an  average  kQ  value  of 
1 . 4  x  10  ^cm  3s  ^ . 
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De-excitation  rate  constant  kg  for  reaction  NO*  +  0 
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APPENDIX  A 

ANALYTICAL  SOLUTIONS  FOR  EXPONENTIAL  n*  OR  n 


The  integral  in  Eq.  (2)  can  be  evaluated  in  terms  of  the  modified 

A 

Bessel  function  of  the  second  kind  when  n  (H)  is  exponential.  One  can  show 
that  the  Bessel  function  reduces  to  an  exponential  if  certain  geometrical 
approximations  are  made.  We  will  derive  this  exponential  solution  starting 
directly  from  the  geometrical  approximations. 


The  partial  derivative  in  Eq.  (2)  is  giver,  by  (see  Fig.  la) 


3z(H,Ht)/3H 


-  (R  +  H)  i(R  +  H)2  -  (K  +  Ht)2] 
a  (R/2)1/2(H  -Hrr1/2 


(A-l) 


where  the  second  form  results  from  expanding  the  terms  in  the  radical  and 
neglecting  H  and  compared  to  R,  the  Earth's  radius.  With  this  approxi¬ 
mation  and  assuming  the  distribution 

n*(H)  =  n"  exp[-(H  -  H0)/H$3  (A-2) 


we  can  write  Eq.  (2)  in  the  form 


n(ht) 


-(ht-h0)/hs  ?  e  -*/Hs 

no e  /  ~rn  dx 


The  definite  integral  is  equal  to  YitH,.;  hence 


(A-3) 


N(HT)  sx  (2A/«*ir)(TrRHs/2)1/2  n*  exp[-(HT-H0)/H$] 
•  (2AA*)URHs/2),/2  n*(HT)  . 


(A-l*) 


The  contribution  of  an  Earth-concentric  layer  (H?,H2)  to  the  radiance 
observed  at  tangent  height  HT  is  obtained  by  generalizing  the  integration 
1 Imits  in  Eq.  (A-3) : 
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N(HT,HrH2)  “  T^^2/  nOe  J  ,  .1/2 


M,  (H"HT) 


dH  .  (a-5) 


This  result  can  be  expressed  In  terms  of  n*(Hy)  and  the  error  function: 


N(Ht,Hj,H2)  -  (2AAir)(irRHs/2),/2n*(HT) 

x  £erf j[(H2-HT)/Hs3,/2|  -  erf  ^ [ ( H 1 -HT) /H$ ] ,/2}J  ; 

H1  >.  Ht  .  (A-6) 


The  limb  profile  given  by  Eq.  (A-1*)  or  (A-6)  must  be  convolved  with 
the  sensor  field  of  view  function  to  obtain  the  observed  limb  radiance 
profile,  unless  n*  varies  by  a  relatively  small  amount  over  the  vertical 
distance  resolved  by  the  sensor. 


APPENDIX  B 

INVERSION  USING  FOURIER  TRANSFORMS 


The  predicted  limb  radiance  profile  for  the  optically -thin  case  is 
obtained  by  evaluating  Eq.  (2)  and  convolving  the  result  N(Hj)  with  the 
sensor  IFOV  function  I (H^) .  If  the  geometric  approximation  represented  by 
Eq.  (A-l)  is  used  in  evaluating  N(Hj),  the  predicted  or  observed  limb 
radiance  profile  can  be  expressed  as 

00 

N(Ht)  at  l(HT)*(2AA*)(R/2)1/2y  dH- (H-HT)‘1/2n*(H)  (B-l) 

HT 

The  integral  in  Eq.  (B-l)  Is  another  convolution;  thus, 

N(Ht)  ~  B-l*W*n*  , 

where 

B  ■  const.  =  (2AMtt)  (R/2)  ^2 

and 

(-H)"1/2  for  hiO 

W  ® 

0  otherwise  . 

uy  the  convolution  theorem,  Eq.  (B-2)  has  the  Fourier  transform 

F(N)  -  BF(l*W) F(n*)  (B-A) 

where  F(  )  derotes  a  Fourier  transform;  likewise,  F  '(  )  will  denote  the 
inverse  Fourier  transform.  The  desired  Inversion  solution  n*(H)  can  be 
obtai ned  as 

n*(H)  F~]  |f(N)/[2F(i*W)]|  .  (B-5) 

The  effective  weighting  function  I*W  is  finite  everywhere  and  can  be 
truncated  (and  apod i zed),  since  n*  •  0  above  some  high  altitude.  Hence, 


(B-2) 


(B-3) 
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A' 

Its  Fourier  transform  exists  and  can  be  evaluated  once  and  for  all. 
Therefore,  the  routine  evaluation  of  Eq.  (B-5)  for  a  given  sensor  involves 
only  two  Fourier  transforms,  which  can  be  performed  using  the  Fast  Fourier 
Transform  (FFT)  method,  and  one  intermediate  division  (multiplication). 

Note  that  F(l*W) ,  and  hence  F (N) ,  can  have  zeroes  at  some  spatial 
frequencies.  This  will  make  n*  the  transform  of  a  function  that  is  indeter¬ 
minate  at  some  frequencies,  but  it  Is  likely  that  the  indeterminancies  can 
be  resolved  by  interpolation. 

We  emphasize  that  this  method  of  limb  radiance  inversion  is  untested 
and  restricted  to  the  optical ly-thin  case. 


*Note  that  F(W)  does  not  exist. 
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